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ABSTRACT 

We present a study of the chemistry of a dense photon-dominated region (PDR) using 
a time-dependent chemical modeL Our major interest is to study the spatial distri- 
bution of complex molecules such as hydrocarbons and cyanopolyynes in the cool 
dense material bordering regions where star formation has taken place. Our standard 
model uses a homogeneous cloud of density 2 x 10^ cm""^ and temperature T = 40 
K, which is irradiated by a far-ultraviolet radiation field of intermediate intensity, 
given by X — 100. We find that over a range of times unsaturated hydrocarbons 
(e.g., C2H, C4H, C3H2) have relatively high fractional abundances in the more exter- 
nal layers of the PDR, whereas their abundances in the innermost layers are several 
orders of magnitudes lower. On the other hand, molecules that are typical of late- 
time chemistry are usually more abundant in the inner parts of the PDR. We also 
present results for models with different density, temperature, intensity of the radia- 
tion field and initial fractional abundances. Our results are compared with both high- 
and moderate-angular resolution observations of the Horsehead nebula. Our standard 
model is partially successful in reproducing the observations. Additional models run 
with different physical parameters are able to reproduce the abundance of many of the 
observed molecules, but we do not find a single model that fits all the observations at 
the same time. We discuss the suitability of a time-dependent model of a dense PDR 
such as ours as an estimator of the age of a PDR, provided that enough observational 
data exist. 

Key words: ISM: molecules - ISM: abundances - Stars: formation - astrochemistry 
- molecular processes 



1 INTRODUCTION 

Far-ultraviolet (FUV) photons (6 eV< hiy < 13.6 eV) can 
drastically affect the structure, chemistry, thermal balance, 
and evolution of portions of the interstellar medium, even 
if only the mean galactic interstellar radiation field (ISFR) 
is present. For this reason, the simplest chemical models of 
dense interstellar clouds should include a treatment of the 
radiative transfer of this radiation through the outer and 
intermediate layers of material, because the field affects ob- 
servations that also sample the interior gas. Indeed, most 
of the molecular gas and certainly all the atomic gas of the 
Galaxy is found in regions where the FUV photon fiux is 
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several times that of the mean ISRF (iHollenbach fc Tielend 
Il997h . Photon-dominated regions (PDRs) are special regions 
of neutral gas exposed to intense fiuxes of FUV radiation, 
which dominate the physical (heating) and chemical pro- 
cesses that take place in their interiors. PDRs can be ob- 
served in very different environments, including dense molec- 
ular clouds such as the Orion Bar or the Horsehead nebula, 
planetary nebulae, winds of red giants and AGB stars, the 
interstellar medium of starburst galaxies such as M82, low 
metallicity dwarf galaxies, etc. PDRs are studied through 
atomic spectral lines such as [C II] at 158/im, [O I] at 63 and 
146 /im, [C I] at 370 and 609 /im; many H2O lines; millimetre 
and sub-millimetre rotational transitions of molecules such 
as CO, ^^CO, or C^^O; ro- vibrational and pure rotational 
lines of H2; and many broad mid-IR features attributed to 
PAHs. 

Physical conditions of PDRs are thus very diverse. They 
can be diffuse, with gas density n '-^ 10 — 100 cm~"^, or 



2 0. Morata & E. Herbst 



Table 1. Standard initial (dark cloud) fractional abundances 
with respect to Jijj 
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dense, with n > 10* cm ^, while the incident FUV flux may 
range from the ISRF to 10^ times the ISRF very near to an 
O star. An important characteristic of PDRs is that they 
have a layered structure, as a result of the interaction of the 
radiation with the gas and dust. Typically, they contain an 
outer layer of partially ionised gas, where H is atomic and 
the predominant form of carbon is CII; a layer of neutral 
gas, where H2 is predominant and carbon is present mainly 
as C I; and a molecular layer, where CO is the main form of 
carbon. The ratio between the flux of the incident radiation 
fleld and the gas density determines the size of these layers. 
The distinct layers should be observable if the PDR is seen 
edge-on, because then the geometry of the cloud is more 
apparent. One of the main difficulties of the study of PDRs 
is, however, the determination of the geometry of the object. 
The gas temperature is also dependent on the strength of 
the radiation field and the position inside the cloud, and can 
range from 10-100 K in the molecular layer to more than 10^ 
K in the surface layers. 

Dense PDRs are of particular interest in our attempt 
to understand the process of star formation. Newly formed 
stars produce high intensity photon fields deep inside clouds 
that lead to the formation of HII regions and, further re- 
moved from the exciting stars, PDRs. The gas-phase chem- 
istry that occurs under these conditions is in most lay- 
ers greatly different from the chemistry found in quiescent 
regions. The abundances of some atoms or molecules are 
greatly enhanced by the presence of the FUV field, while 
other species are destroyed by the radiation. For instance, 
the abundances of several hydrocarbons (C2H, C-C3H2, C4H, 
CeH) have been found to be almost a s high in relatively 
unshielded regions as in dark clouds (|Fuente et al.l l2003l : 
iTevssier et ai] |2004| ). However, some of these high abun- 
dances are difficult to explain with normal gas-phase chemi- 
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Figure 1. Fractional abundances of H (filled triangles), H2 (filled 
squares), C (filled circles), C"*" (filled triangles pointed down- 
ward), and CO (filled diamonds) with respect to nn as a function 
of Ay for our standard model. Three different times are shown: 
10* yr (dashed-dotted lines), 4 X 10* yr (dashed lines), and 5 X 10^ 
yr (solid lines). The solid lines correspond to steady-state abun- 
dances. 



cal models and a lin k between small hydr ocarbons and PAHs 
has been proposed (|Tevssier et al.|[2004l ). 

In this paper, we study the chemistry of a dense PDR 
created after a star has been formed inside a dark cloud. 
We are particularly interested in the spatial distributions 
of hydrocarbons and other complex molecules, and whether 
or not time-dependent chemistry can explain some of the 
discrepancies found between the obse rved abundanc e s and 
the results of steady-state models (|Tevssier et al. I l2004l : 
iPetv et al.ll2005l ) for the Horsehead nebula, as well as pro- 
viding us with some estimate of the age of this PDR. An- 
other use o f time -dependent PDR chemistry is described by 
iBell et al.l (|2005| ). The structure of the remainder of the 
paper is as follows: in Section [2] we describe the time- 
dependent chemical model used. Our results obtained with 
our standard model are then considered in Section |31 while 
in Section |4] we discuss the dependence of our results on as- 
sorted physical parameters. The agreement between model 
results and observations of the Horsehead nebula is probed 
as a function of time in Section[5] Finally, Section[6]contains 
a summary of the paper. 



2 MODEL 



In our model, based on lLee et all (|l996l ). we assume a semi- 
infinite plane-parallel cloud divided into rii slabs, on one side 
of which impinges an ultraviolet radiation field, defined in 
our case in units of the iDraind l| 19781 ) standard radiation 
field scaled by a multiplicative factor x- The visual extinc- 
tion, Av, and the total col umn density, Nh, are related b y 
Kv = 6.289 X 10"^^ Nh (jWagenblast fc Hartguistl Il989l ). 
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Figure 2. Fractional abundances of the molecules CH4 (filled 
diamonds), C2H (filled circles), C4H (filled squares), and C3H2 
(filled triangles) with respect to njj as a function of Ay for our 
standard model. Three different times are shown: IC* yr (dashed- 
dotted line), 3 X 10^ yr (dashed line), and 5 X 10^ yr (solid line). 
The solid lines correspond to steady-state abundances. 



The model treats radiative transfer and, in particular, takes 
into accou nt the self- and cross-shield ins due to H2 and CO. 
We use the lDraine fc Bertoldil (|l996l ) analytical approxima- 
tion formula to calculate the H2 photo-dissociation rate (fe 
= 1 km s~^), which we have found to be in closer agreement 
with the results of the 'exact' treatm ent of H2 self-shielding 
llLe Petit. Roueff fc Le Bourlotll200^ than the method pre- 
sented in Lee et al. ( 19961 ). However, the larger differences 
are only found for small values of the visual extinction, 
Ay < 0.01, which are not the main source of interest of 
the present paper. For th e CO photodisso ciation rate, we 
use the same method as in iLee et al.1 (Il996l ). 

For our standard model, we have assumed a homoge- 
neous cloud with physical conditions similar to those ob- 
served in the Horsehead nebula (see Section [5}, namely 



a density nu ~ 2 x 10 cm and a temperature T — 
40 K. The cosmic-ray ionisation rate, ^, is assumed to be 
1.3 X 10"" s-\ We use a value for the scaling factor of 
the radiation field of x = 100, which is similar to the 
one measured in th e Horsehead nebula (|Zhou et al.l 1 19931 : 
lAbergel et al]|2003l ). In our standard model, the cloud is 
divided into 42 slabs of differing length, defined in a way 
so as to have a smooth distribution of visual extinction 
points from Ay = 5 x 10~* to Ay = 13, but good cov- 
erage of the positions where the fractional abundances of 
most of the species show steeper slopes. The time-dependent 
gas-phase chemistry is solved slab-wise, starting from the 
outermost slab and moving successively inward. The model 
solves a system of ordinary differential equations with the 
Gear method for each slab using the H2 and CO photodis- 
sociation rates from the (i-l)th slab as the input values 
for the i-th slab. For the che mistry, we use the osu.2003 
jSmith. Herbst fc Chan3l2004l ) gas-phase chemical network. 




Figure 3. Same as Fig. [21 but for the molecules C3H (filled 
squares), CgH (filled circles), and C3H4 (filled triangles). 



which contains 419 species and 4103 reactions. This rela- 
tively old version of the osu network was used for a variety 
of reasons. First, it was the model with which the origi- 
nal calculations were begun. During the calculations, sev- 
eral newer versions of the network (osu. 2005) were pre- 
pared, but required modifications and corrections, which 
were not fully resolved until the summer of 2007. To de- 
termine if the newest version of the osu code (osu_03_2008; 
see http:/ /www. physics. ohio-state.edu/'^eric/ ) leads to any 
significant changes, we have run a new version of our stan- 
dard model with this network. The differences are discussed 
in Appendix A. 

We assume that our dense PDR is formed as a result 
of the process of star formation inside a dense cloud, when 
radiation from a new star begins to affect the chemistry of 
the surrounding gas. Thus, we assume that, at t = 0, the 
gas has already undergone chemical evolution during ~ 10^ 
yr under dense conditions: a density nn = 2 x lO'' cm~'^ 
and T = 10 K. The initial abundances are those of a typical 
dark cloud (see Table [l| as derived from the values measured 
for s maller species (comp osed of up to five atoms) in TMC- 
ICP (jSmith et al.l|2004| ). For those smaller species where 
there is no available observed v alue, we ado p ted t he frac- 
tional abundances calculated bv lSmith et al. I (|2004l ) at lO'^ 
yr. In our model, we maintain the standard depleted cosmic 
abundance ratios C/H= 7.3 x lO"'' and 0/H= 1.76 x 10"*. 
Since the total C/H abundance measured in TMC-IP is a 
factor ^ 2 less than the standard cosmic ratio, we have as- 
sumed that the remaining C is locked up in CO, with the 
initial abundance shown in Table [T] The initial abundance 
of atomic O is then chosen so that the overall standard de- 
pleted O /H abundance is used. We discuss in Section 14.31 
other possible results using the lower initial abundances of 
CO. 

The t = Q values listed in Table [1] are more general 
than those found for TMC-1, where a much more complex 
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Table 2. Peak fractional abundances and for selected hydrocarbons in the standard model 
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molecule development has already occurred, which appears 
to be unique. 



3 RESULTS 

Fig. [1] shows the fractional abundance profiles of H, H2, C, 
C"*", and CO obtained with our standard model at selected 
times up to steady state (5 x 10^ yr). Such a profile at t = 
would be totally fiat at the value given in Table 1. Here the 
three layers described in Sect. [2] can be distinguished. The 
crossing point between the fractional abundances of H and 
H2 moves inwards with time, and it is found at Av ^ 0.04 
at steady-state. The low value for this visual extinction is 
somewhat artificial; proper inclusion of factors such as a 
lower density and higher temperature at the cloud edge and 
a less efficient surface formation of H2 would increase the ex- 
tinction at crossing. Further into the cloud, CO becomes the 
most abundant carbon-bearing species at Ay ~ 1.7, while 
C becomes more abundant than C"*" at visual extinctions 
between 3 and 8.5. 

3.1 Hydrocarbons 

Figs [2] and [3] show the fractional abundances of several hy- 
drocarbon molecules as a function of Av for our standard 
model at three different times: t = 10*, 3 x 10^, and 5 x 10® 
yr. The first two times correspond, respectively, to an early 
time in the evolution when only the outermost layers are 
approaching steady state, and a time when steady-state is 
virtually reached in the outer to intermediate regions of the 
FDR {Av < 3 — 5). In addition to the figures, Table[2]shows 



the peak fractional abundances and corresponding visual ex- 
tinctions of most hydrocarbons present in our model at the 
three different times. In this and subsequent tables of frac- 
tional abundances, the order of presentation is in terms of 
decreasing steady-state peak fractional abundance. 

The evolution of the fractional abundances of a hydro- 
carbon depends upon whether the molecule is present ini- 
tially in the gas, a group that includes C2H, C4H, and C3II2, 
or whether it is not, a group that includes both complex 
species and hydrogenated species such as CH4 and C3II4. 
As emphasized in Fig. 3] the main difference is that the 
molecules pertaining to the first group (C2II in Fig. |4]) go 
through an early phase of selective destruction before reach- 
ing steady-state, while the molecules belonging to the second 
group (CH4 in Fig. |4]) start with a phase of formation. 

Looked at in more detail, the unsaturated hydrocarbons 
initially present in the gas are quickly heavily destroyed in 
the very outer layers of the FDR, while the destruction is 
much slower in more shielded layers. By a time of a fewx 10* 
yr, fractional abundance peaks are found at Av ~ 1.5-2 
with a value close to that of steady-state. As time pro- 
gresses, steady-state is reached for the outer parts, while 
the fractional abundances in the inner layers may drop sev- 
eral orders of magnitude before steady-state is reached at 
t ~ 5-10 X 10® yr. The end resuh is a large decrease in 
abundance from the initial value for the outermost layers 
(Ay < 1) and a significant decrease, very much dependent 
on the molecule, for slabs with Av > 7. 

The molecules that were not present initially in the gas 
(e.g. CH4, CII2) initially show a larger increase in the deeper 
layers of the cloud. We usually find the highest abundance 
for each Av at a time t ~ 10* yr or shortly thereafter. After 
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Figure 4. Fractional abundances with respect to rijj for the 
molecules C2H (left panel) and CH4 (right panel) as a function 
of Av in our standard model for five different times: 10^, 10*, 
10^, 10^, and 10^ yr. The solid lines represent steady-state abun- 
dances. 



this moment, the evolution of the fractional abundances of 
these molecules is similar to the one described for the previ- 
ous group: the fractional abundances do not change appre- 
ciably in the outer parts of the PDR (Av ~ 1-2) while the 
abundances may drop, although not so dramatically, in the 
inner layers until reaching steady-state conditions at t ~5" 
10x10® yr. Interestingly, CH4 abundance does not follow 
this description completely for Av > 3, where it keeps in- 
creasing until a time 5 x lO'* yr. Then, it shows very little 
variation, and it only decreases very slightly for Av > 8. 

All hydrocarbons show a similar abundance profile at 
steady state (or even at earlier times): their fractional abun- 
dances grow with depth into the cloud until they reach a 
peak, after which the fractional abundance decreases with 
depth more or less drastically. The peak position and the 
steepness of the decrease depend on the hydrocarbon. As can 
be seen in Table [2] the more unsaturated and simpler hy- 
drocarbons typically peak in outer layers of the cloud, while 
as the number of C atoms grows or as the molecules become 
more hydrogenated, the peaks are found deeper into the 
cloud. For instance, by steady state, C2H peaks at Av ~ 1.6; 
C3H2, at Av = 1.9; CgH, at Av = 2.0; and CH4, at Av = 8. 
This distribution is not continuous; there seem to be two 
main regions, Av = 1.5-2 and Av = 5.5-8, where all the 
peaks are found. In fact, most of the hydrocarbons show two 
peaks at those two ranges, the relative importance of which 
seems to depend on the number of carbon atoms present or 
on being or not saturated. This secondary peak (primary in 
some cases) is due to the bump in C abundance shown in 
Figure 1. Hydrocarbons containing fewer carbon atoms typ- 
ically reach larger fractional abundances; for example, CH4, 
C2H and C4H are some of the most abundant hydrocarbons. 



3.2 HCN, HNC, and cyanopolyynes 

Fig. [5] shows the evolution of the fractional abundance 
profiles for the HCN and HNC molecules, and the 
cyanopolyynes HC3N, HC5N, HC7N, and HCgN, while Ta- 
ble [3] lists the abundance peaks and the Av where they are 



>< 




Figure 5. As Fig. [J] but for the molecules HCN (diamonds), 
HNC (stars), HC3N (filled triangles pointing down), HC5N (filled 
squares), HC7N (filled triangles pointing up), and HCgN (filled 
circles). 



found. Initially present in the gas, HC3N experiences a pro- 
gressive diminution in all the slabs as time progresses, which 
is particularly severe in the first 10^ yr. When steady-state 
is reached, HC3N shows a profile with two peaks, a very low 
one at Av ~ 2, and the main one at Av ~ 8. Unlike the 
case for hydrocarbons, the abundance of HC3N in the deep- 
est layers is less than an order of magnitude below that at 
the peak. 

The more complex cyanopolyynes show a behaviour 
similar to those unsaturated hydrocarbons not present ini- 
tially in the gas. First, they go through a phase of forma- 
tion until they achieve maximum abundances at all Av at 
a time t ^ 10* yr, when peak abundances are in the inner 
regions. Then they are progressively destroyed at all optical 
depths. The steady-state profiles of these molecules show a 
low abundance peak at Av ~ 2, the main abundance peak, 
one or two orders of magnitude over the secondary peak, at 
Av ~ 4 — 8, and a decrease at deeper layers. 

The initially-present HCN and HNC molecules show a 
completely different but nearly identical behaviour (see Ta- 
ble[3]|. Their destruction is very intense in the first 10* yr for 
slabs Av < 2, while only moderate for the deeper layers. The 
destruction goes on until a time t 10^ yr, when there is a 
slight increase of abundances for Av ~ 3-8, before steady- 
state is reached. An abundance peak then exists at Av ~ 5.5 
and a slightly lower, but almost uniform, abundance occurs 
at deeper layers. As happened with the cyanopolyynes, there 
is a small bump in abundance at Av ~ 2.1, approximately 
20 times lower than the main peak. 

3.3 Other molecules 

Fig-Elshows the abundance profiles at assorted times of CS, 
HCO^, NH3, and N2H^, all of which are well-known species 
in dense regions. The peak abundances at the selected times 
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Table 3. Peak fractional abundances and Ay for cyanopolyynes, HCN, and HNC 



Species 


t = 10* 




t = 3 X 10^ 




t = 5 X 10*^ 


yr 




n(X)/nH 


Av 


n{X)/nH 


Av 


n(X)/nH 


Av 


HCN 


3.9 X 10"'' 


13.0 


1.0 X 10"^ 


8.0 


1.5 X 10-» 


8.0 


HNC 


3.5 X 10-s 


13.0 


9.9 X 10-"' 


9.0 


1.4 X 10-s 


8.0 


HC3N 


4.7 X 10"" 


13.0 


3.3 X 10-1^ 


8.0 


2.4 X 10-1^ 


8.0 


HC5N 


1.1 X IQ-^^ 


13.0 


4.5 X 10-" 


13.0 


5.4 X 10-15 


6.0 


HC7N 


4.5 X 10-1* 


13.0 


2.4 X 10-15 


13.0 


5.3 X 10-1*5 


5.5 


HC9N 


1.7 X 10-15 


13.0 


1.8 X 10-1"^ 


6.0 


6.7 X 10-1^ 


5.5 



Table 4. Peak fractional abundances and Av for other molecules typically detected in dense clouds 



Species 


t = lO"! 


yr 


t = 3 X 105 


yr 


t = 5 X 10*5 


yr 




n{X)/nH 


Av 


n{X)/nH 


Av 


n{X)/nH 


Av 


NH3 


6.8 X lO-'J 


13.0 


1.7 X 10-8 


13.0 


3.6 X 10-8 


13.0 


HCO+ 


4.5 X 10"'' 


13.0 


8.3 X 10-9 


13.0 


9.4 X 10-9 


13.0 


H2CO 


1.2 X 10-9 


13.0 


1.1 X 10-9 


7.0 


1.1 X 10-9 


7.0 


CS 


1.5 X 10-'-* 


13.0 


6.3 X 10-10 


6.0 


5.9 X 10-10 


6.0 


N2H+ 


1.0 X 10-1" 


13.0 


1.9 X 10-10 


13.0 


3.4 X 10-10 


13.0 


HCS+ 


1.3 X 10-11 


2.0 


1.5 X 10-11 


2.2 


1.5 X 10"ii 


2.2 




Figure 6. As Fig. [2] but for the molecules CS (filled squares), 
HCO+ (filled diamonds), NH3 (filled circles), N2H+ (filled trian- 
gles pointing up), and HCS+ (filled triangles pointing down). 



are listed in Table |4] Present in the initial gas, CS shows 
a fractional abundance profile quite similar to unsaturated 
hydrocarbons with two peaks of very similar abundances, 
the inner one ~ 1.5 times larger than the outer one, but 
with a very small fall in abundance for the deeper layers, 
where abundance is very close to constant. 

The species HCO'*', NH3, and N2H^ show profiles some- 



what similar to those of HCN and HNC. Unlike these latter 
two species, however, the fractional abundance peak of the 
three molecules is found in the innermost slab of our model. 

One common thread running through these descrip- 
tions of assorted profiles is the following: the dominant fea- 
tures of the fractional abundance profiles as a function of 
Av for most species show a strong similarity to plots of 
the fractional abundance of the same molecules as a func- 
tion of time in hom ogeneous dark cloud models (see e.g., 
iTavlor. Morata fc WiUiama IQga ). The so-called 'early-time' 
molecules, such as CS or C2H, peak in the outer layers of the 
PDR, while the 'late-time' molecules, such as NH3, N2H+, 
or HCO"*", peak in the innermost layers. This is not too sur- 
prising; the early-time chemistry in the homogeneous cloud 
(based on a high abundance of atomic C) is similar to that 
in the PDR layers with Av ~l-2, while the late-time chem- 
istry (based on a high abundance of CO) is similar to that 
in the inner layers. 

4 INFLUENCE OF PHYSICAL PARAMETERS 
4.1 Radiation field 

In order to determine how the external radiation field af- 
fects the abundances of the different species, especially the 
hydrocarbons, we ran three additional models, each with a 
different intensity of the radiation field, indicated by scaling 
factors of X = 10, 10^, and 10*. Fig.[7]compares the results of 
our standard model with the results obtained with a model 
using X = 10* at an early stage of evolution, 10* yr, and 
at 10^ yr, when steady-state has been reached for all Av. 
Table [5] lists the peak abundances of a variety of species for 
three times when x = 10*. 
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Table 5. Peak fractional abundances and for selected molecules for a model with x = 10* 



Species t = lO" yr t = 3 X 10^ yr t = lx lO'^ yr 

n(X)/nH Av^ n{X)/nH A^ n(X)/nH Ay 



CH4 


5, 


.7 


X 


10" 


-9 


9, 


,0 


2, 


,3 


X 


10- 


-8 


13, 


,0 


2, 


,3 


X 


10- 


-8 


13, 


,0 


NH3 


6, 


.7 


X 


10- 


-9 


13, 


,0 


1, 


1 


X 


10- 


-8 


13, 


,0 


1, 


,9 


X 


10- 


-8 


13, 


,0 


CH 


6, 


.0 


X 


10- 


-9 


2, 


,5 


1, 


,2 


X 


10- 


-8 


3, 


,0 


1, 


,2 


X 


10- 


-8 


3, 


,2 


HCO+ 


4, 


.4 


X 


10- 


-9 


13, 


,0 


6, 


,8 


X 


10- 


-9 


13, 


,0 


7, 


1 


X 


10- 


-9 


13, 


,0 


CH2 


4, 


.6 


X 


10- 


-9 


2, 


,5 


5, 


,4 


X 


10- 


-9 


3, 


,2 


5, 


,5 


X 


10- 


-9 


3, 


,2 


C2H 


8, 


.1 


X 


10- 




2, 


,5 


3, 


,3 


X 


10- 


-9 


3, 


,0 


3, 


,3 


X 


10- 


-9 


3, 


,2 


C2H2 


4, 


,1 


X 


10- 


-9 


13, 


,0 


1, 


,7 


X 


10- 


-9 


13, 


,0 


1, 


,6 


X 


10- 


-9 


13, 


,0 


HCN 


3, 


.8 


X 


10- 


-9 


13, 


,0 


1, 


,1 


X 


10- 


-9 


13, 


,0 


1, 


,6 


X 


10- 


-9 


13, 


,0 


HNC 


3, 


.5 


X 


10- 


-9 


13, 


,0 


1, 


,0 


X 


10- 


-9 


13, 


,0 


1, 


,5 


X 


10- 


-9 


13, 


,0 


H2CO 


1, 


,2 


X 


10- 




13, 


,0 


1, 


1 


X 


10- 




9, 


,0 


1, 


1 


X 


10- 




9, 


,0 


CS 


1, 


,5 


X 


10- 


-9 


13, 


,0 


9, 


,0 


X 


10- 


-10 


9, 


,0 


8, 


,7 


X 


10- 


-10 


9, 


,0 


N2H+ 


9, 


.9 


X 


10- 


-11 


13, 


,0 


1, 


,6 


X 


10- 


-10 


13, 


,0 


2, 


,6 


X 


10- 


-10 


13, 


,0 


CH3 


9, 


.2 


X 


10- 




7, 


,0 


1, 


,2 


X 


10- 


-10 


6, 


,0 


1, 


,2 


X 


10- 


-10 


6, 


,0 


C3H2 


1, 


.0 


X 


10- 


-9 


13, 


,0 


7, 


1 


X 


10- 


-11 


9, 


,0 


7, 


,3 


X 


10- 


-11 


9, 


,0 


C3H 


1, 


.9 


X 


10- 


-10 


13, 


,0 


5, 


,5 


X 


10- 


-11 


9, 


,0 


5, 


,6 


X 


10- 


-11 


9, 


,0 


C4H 


8, 


.7 


X 


10- 




13, 


,0 


1, 


,6 


X 


10- 




9, 


,0 


1, 


,6 


X 


10- 




9, 


,0 


C2H4 


3, 


.9 


X 


10- 


-13 


13, 


,0 


2, 


,1 


X 


10- 


-12 


13, 


,0 


2, 


,4 


X 


10- 


-12 


13, 


,0 


C2H3 


1, 


.4 


X 


10- 


-12 


7, 


,0 


2, 


,2 


X 


10- 


-12 


8, 


,0 


2, 


,2 


X 


10- 


-12 


8, 


,0 


C4H2 


2, 


,3 


X 


10- 


-10 


13, 


,0 


2, 


,5 


X 


10- 


-12 


13, 


,0 


1, 


,9 


X 


10- 


-12 


13, 


,0 


CsH 


1, 


.0 


X 


10- 


-11 


13, 


,0 


8, 


,3 


X 


10- 


-13 


9, 


,0 


7, 


,6 


X 


10- 


-13 


9, 


,0 


CeH 


4, 


.1 


X 


10- 


-12 


13, 


,0 


7, 


,9 


X 


10- 


-13 


9, 


,0 


7, 


,0 


X 


10- 


-13 


9, 


,0 


HC3N 


4, 


,5 


X 


10- 


-11 


13, 


,0 


3, 


,4 


X 


10- 


-13 


13, 


,0 


2, 


,5 


X 


10- 


-13 


13, 


,0 


C3H4 


9, 


,8 


X 


10- 


-14 


13, 


,0 


2, 


,3 


X 


10- 


-13 


13, 


,0 


2, 


,0 


X 


10- 


-13 


13, 


,0 


CgH 


2, 


,0 


X 


10- 


-13 


13, 


,0 


1, 


,3 


X 


10- 


-13 


9, 


,0 


1, 


1 


X 


10- 


-13 


9, 


,0 


C7H 


4, 


.4 


X 


10- 


-13 


13, 


,0 


1, 


1 


X 


10- 


-13 


9, 


,0 


8, 


,5 


X 


10- 


-14 


9, 


,0 


C3H3 


1, 


.3 


X 


10- 


-14 


13, 


,0 


3, 


,4 


X 


10- 


-14 


9, 


,0 


4, 


,4 


X 


10- 


-14 


9, 


,0 


C5H2 


5, 


.1 


X 


10- 


-12 


13, 


,0 


2, 


,9 


X 


10- 


-14 


13, 


,0 


2, 


,8 


X 


10- 


-14 


9, 


,0 


C4H3 


1, 


.3 


X 


10- 


-13 


13, 


,0 


2, 


1 


X 


10- 


-14 


9, 


,0 


2, 


,3 


X 


10- 


-14 


9, 


,0 


CH3C4H 


1, 


.1 


X 


10- 


-13 


13, 


,0 


1, 


,9 


X 


10- 


-14 


13, 


,0 


1, 


,4 


X 


10- 


-14 


9, 


,0 


HC5N 


1, 


,1 


X 


10- 


-12 


13, 


,0 


2, 


,6 


X 


10- 


-14 


13, 


,0 


1, 


,3 


X 


10- 


-14 


9, 


,0 


C6H2 


5, 


,3 


X 


10- 


-12 


13, 


,0 


1, 


,9 


X 


10- 


-14 


13, 


,0 


1, 


,2 


X 


10- 


-14 


9, 


,0 




4, 


,7 


X 


10- 


-14 


13, 


,0 


1, 


1 


X 


10- 


-14 


9, 


,0 


3, 


,2 


X 


10- 


-15 


9, 


,0 



The evolution of the fractional abundances when x = 
10* differs most strongly from the standard the 
right-hand side of Fig. [7] shows, in the sweeping of most 
of the molecular species from the outer layers by the intense 
radiation. Most of the molecules shown have already disap- 
peared from the layers with Ay < 1 hy t — 10'* yr and, 
when steady-state is reached at 10^ yr, there are no sig- 
nificant molecular abundances at slabs with Ay < 1. Even 
molecular hydrogen disappears from the external layers as 
the position at which H2 becomes the most abundant form 
of hydrogen increases to ~ 1.3 when x = lO"- Thus, for 
the model with x = 10 , the gas at Av < 1 is mostly in 
atomic form. The fractional abundance peaks of the species 
with profiles that peak in the outer layers are consequently 
found deeper into the cloud at all times, the typical effect 
being an increase from Ay 1.5-2 to Av ~ 3-4 (and the 
inner maxima and/or enhancements from A^ ~ 4.5-7 to 
Ay ~ 7-9). On the other hand, the results for the model 
with a less intense radiation field than the standard model, 
X = 10, are the opposite and many molecular species have 
significant fractional abundances at Ay as low as 0.1. The 
fractional abundance peaks are then located in more exter- 
nal layers of the cloud. 



The maximum fractional abundances are also modified 
by the intensity of the radiation field in all of the cases. As 
X is increased from 100 to 10*, the outer peak abundance 
of hydrocarbons is typically reduced by a factor of < 10, 
but the decrease can range from very little, as for CH2, to 
two (C4H or CeH) or three (C8H) orders of magnitude. On 
the other hand, the 'inner' abundance peak of hydrocarbons 
increases by a small factor, between 2 and 10, and in some 
cases, such as C3H2, it can become larger than the (outer) 
peak abundance in the standard model. 

The peak fractional abundances of the cyanopolyynes 
have a much greater dependence on the radiation field, and 
we find general enhancements of two and three orders of 
magnitude for the outer layers, when x = 10*. This increase 
does not change the fact that there are only traces of these 
molecules in these layers. The peak abundance at high visual 
extinction also increases with a larger radiation field, but 
only by small factors (2-3). 

The radiation field is able to modify the chemistry even 
in deep slabs of the cloud. We find that even if the changes 
in intermediate layers {Ay ~ 5-8) can be small, a factor 
of a few, they are non-zero; only the gas at Av ^ 9-10 is 
completely shielded from the radiation for all the models 
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Figure 7. Comparison of the fractional abundances obtained 
with (left-hand panel) the standard model and (right-hand panel) 
with a model using x = 

Iff*. The results are shown at an early 
time, 10* yr (dashed lines and open symbols), and at lO'^ yr (solid 
lines and filled symbols), when steady-state has been reached. 



and does not show any significant changes at any time be- 
tween the different models. As a result, for the two radiation 
fields considered, all the molecules that reach their maxi- 
mum abundances in the inner cloud, such as HCN, NH3, and 
N2H''', have approximately the same abundances at high vi- 
sual extinction at steady-state and at intermediate times. 
On the other hand, these same molecules, which were not 
abundant in the outer layers of the cloud for x = 100, com- 
pletely disappear from the superficial layers for x = 10*. 

4.2 Density and temperature 

In order to explore the effects of different densities and tem- 
peratures, we reran our standard model, modifying either 
the density or the temperature of the gas. The densities and 
temperatures used were nn ~ , 10^, and 10^ cm"'', and 
T = 10, 100, and 200 K. 

The effects of density change bear some similarities with 
the effects of the radiation field: a higher density increases 
the shielding and the abundances of most of the molecules 
in the outer layers of the cloud; a lower density on the other 
hand, clears these same layers of much molecular material, 
and displaces the peak abundance position to deeper parts of 
the cloud. The results are not identical, particularly in the 
way the inner and outer peaks of many molecules change 
with density. 



Table 6. Atomic initial abundances 



Element n(X)/nH 



H 


1 


He 


0.14 


C+ 


7.3 X 10-5 


N 


2.14 X 10-5 





1.76 X 10-* 


S+ 


8 X 10-* 


Si+ 


8 X 10-^ 


Fe+ 


3 X 10-8 


Na+ 


2 X 10-8 


Mg+ 


7 X 10-8 


P+ 


3 X 10-8 


C1+ 


4 X 10-8 



Unsaturated hydrocarbons (C2H, C4H, CeH, CgH) 
show an increase in the peak abundance in the outer abun- 
dance peak with increasing density, that can range from a 
few percent (C2H) to an order of magnitude (C4H, CeH); 
and a simultaneous decrease in the abundance of the in- 
ner abundance peak, which can be almost three orders of 
magnitude for the molecules containing more than a few 
carbon atoms. The resulting effect is that at high densities 
there is only a narrow range of Ay in the outer layers of 
the cloud where these molecules show a significant abun- 
dance; the more carbon atoms they contain, the narrower 
the range. At the lowest density, on the other hand, the in- 
ner abundance peak may not only be larger than the outer 
abundance peak, but also larger than the peak abundance 
in the standard model. Hydrocarbons that peak inside the 
cloud, such as methane and C3H4, also show a clear in- 
crease in abundance as density increases. At n = 10^ cm-'', 
methane also shows an outer peak. 

The rest of the molecules we have discussed in Section [3] 
also show clear trends in peak abundance with a change 
in density. The cyanopolyynes show a steep increase in the 
abundance in the outer layers as density increases (it can 
reach three or more orders of magnitude) , but never exceed 
the inner peak abundance, which does not change much and 
is found at slightly lower Ay. At the higher density, HC3N 
reaches lO-*'^ cm-'' for times ~ 5x 10^ yr. The abundance of 
the radical CS increases slightly with density (up to a factor 
of 4 in the outer parts of the cloud), whereas HCO'*' shows 
a behaviour similar to the unsaturated hydrocarbons, but 
with changes typically less than an order of magnitude: the 
peak abundance at high density is a factor ~ 8 lower than 
for the standard model. NH3 and N2H"'" show similar trends 
as HCO"'', but while the fractional abundance decrease at 
n = 10^ cm-^ is less than a factor of 3 for NH3, it is larger 
than 30 for N2H+. 

Variations in the peak fractional abundance with tem- 
perature as we raise T to 200 K tend to be less dramatic 
than with density. The differences in fractional abundance 
in the outer layers relative to our standard model are rarely 
larger than a factor of 2, and often show less than a 50 per 
cent change. The differences are larger for the inner layers, 
where a few molecules show a change of more than an or- 
der of magnitude for T — 200 K. Most of the carbon chains 
(except C2H or CeH) show a large decrease of abundance 
(~ 100) in the innermost layers as T increases, whereas for 
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Figure 8. As Fig.[T]for a model with atomic initial abundances. 



the outer layers only C2H, CsH, and CsH show an increase in 
abundance with temperature. The fractional abundances of 
CH4, H2CO, CS, and the cyanopolyynes HC3N and HC5N 
also diminish with increasing T, while the abundances of 
HCN and HNC increase in the outer layers with T (a factor 
of ~ 3 at r = 200 K), but decrease in the inner layers (less 
than 40 per cent). The ions HCO^ and N2H^ show increas- 
ing abundance with increasing temperature (a factor of ~ 2 
for T — 200 K) . The abundances of the rest of the molecules 
that we discussed in Section |3] show small changes of a few 
percent. 




4.3 Initial abundances 

The previous models were all run assuming as initial abun- 
dances those expected if a dense PDR is created by a star 
recently formed inside a molecular cloud. In order to ex- 
plore the dependence of the model on the initial fractional 
abundances, it is also interesting to determine what are the 
differences between our standard model and other models 
that have different initial fractional abundances. 

First, we compare our standard model with a model 
starting from gas that is atomic in nature. Figures [8] 
and [5] show some results obtained with a model identical 
to our standard model, but using the initial abundances 
listed in Table 6l de rived from 'low metal' abundances 
llRufBe fc HerbstBoOOl '). The major difference between Fig.\E\ 
and its analogous. Fig. [TJ is that, with standard initial 
abundances, H2 decreases at low Av as a function of time 
whereas, with atomic initial abundances, H2 must be syn- 
thesized and increases at all Ay as a function of time. Fig. [9] 
clearly shows that the different initial abundances determine 
the early evolution of the gas. The first years are obviously 
spent in the process of building more complex molecules 
from the available elements, while, in our standard model, 
a general destruction of initially present molecules is pre- 
dominant for very early times. Table [7] lists the maximum 
abundances and peak positions of several interesting species 



Figure 9. Fractional abundance profiles for a model with atomic 
initial abundances and x = 100 at three times: 10* yr (dot-dashed 
lines), 3 X 10^ yr (dashed-lines) and 5 X 10^ yr (solid lines). Some 
hydrocarbons are shown on the top figure; viz., C2H (filled cir- 
cles), C3H2 (filled triangles up), and CH4 (filled diamonds). On 
the bottom figure are depicted HC3N (filled circles), NH3 (filled 
triangles up), and CS (filled diamonds). 



at selected times. Much lower than in the standard model 
for a time t = 10'' yr, the fractional abundances show an 
increase as Av grows until a plateau is reached. By 1-3x10^ 
yr, the fractional abundances of some molecules, such as 
many carbon chains and cyanopolyynes, can be much larger 
at Av > 3 than for the standard model at the same times, 
from factors of 2-3 to more than ~ 1000. Eventually, after 
a few xlO^ yr when most of the C begins to be locked in 
the CO molecule, the abundances for many molecules in the 
inner layers follow the same evolutionary path as in the stan- 
dard model and start decreasing, while there is still a slight 
increase in the molecular abundances in the outer layers. The 
net effect is a change from the 'plateau' profile to the same 
peaked profiles we observe for most species in our standard 
model. Interestingly, the abundances of HCO'*', NH3, and 
N2H''', which peak in the innermost layers, are smaller than 
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Table 7. Peak fractional abundances and for selected molecules for a model with atomic initial abundances and x = 100 



Species t = 10"* yr t = 3 X lO''^ yr t = 5 X lO"^ yr 

n{X)/nH Av n(X)/nH Ay n{X)/nH Ay 



NH3 


4, 





X 


10" 


-15 


13 


.0 


7. 


7 


X 


10" 


-9 


13. 





3. 


,6 


X 


10" 


-8 


13. 


.0 


CH4 


1. 


1 


X 


10" 
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the ones in the standard model for times t < 1 — 2 x 10 
yrs. Steady state is reached at f ~ 10^ yr at all Ay- The 
abundances at steady state are the same or almost identical 
(< 1 per cent) to the fractional abundances of the standard 
model. 

As we pointed out in Section O the CO initial frac- 
tional abundance of our standard model is greater than the 
one measured in TMC-IP by a factor ~ 2. In order to see 
how this can affect our conclusions, we ran t wo additional 
models that used the CO abundance listed by ISmith et alj 
("2004') for TMC-1: i) a model (Model CO-1) maintaining 
the total standard elemental amount of C and O, and ii) a 
model (Model CO-2) not maintaining either the C/H ratio of 
7.3 X lO"'^ or the 0/H ratio of = 1.76 x 10"*. Table[8]shows 
the initial abundances of C, O, and CO for these models, as 
well as the overall C/H, O/H, and C/0 elemental abundance 
ratios, designated with the subscript 'e/em'. As an example. 
Fig. [To] shows the C2H abundances in these two models rel- 
ative to the C2H abundance in our standard model. 

For the first model, the more abundant initial free C 
in the gas produces a smaller abundance of CO for times 
< 2 X 10^ yrs and a correspondingly larger abundance of C 
and C^ up to 3-5x10^ yrs. This effect is particularly im- 
portant for atomic C in the inner layers. For early times, 
the free carbon produces an increase of the abundance of 
all the carbon-chain molecules, at all layers Ay ^ 2. These 
differences in abundance relative to the standard model di- 
minish with time from the outer to the inner layers, and last 
up to 5-10x10^ yrs for layers Ay >3 — 4. The differences 



Table 8. Initial abundances of C, O, and CO [n(X) /nji] for var- 
ious models 



species Standard Model Model CO-1 Model CO-2 
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-4 


1.43 X 10" 
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0.41 




0.41 




0.28 





in abundance at the positions of the external peak abun- 
dances of the standard model are generally very small. For 
early times, the abundances at layers with 2<Ay <5 can 
be significantly larger by factors between a few (for CH4) 
and up to 10^ — 10* (for CgH), while for layers with Ay > 5, 
most of the carbon-chain molecules show abundances that 
range from ~ 25 to a few hundred times larger than the 
ones in the standard model at early times. Other carbon- 
containing molecules show a similar behaviour. On the other 
hand, HCO"*" and NH3 show smaller abundances, by factors 
of 3-4, for layers Ay < 6. For times later than ~ 10® yr, 
the abundances are practically identical to the ones from 
the standard model, as we expected from the results of the 
model with initial atomic abundances. 

Although the late-time abundances for the second 
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Figure 10. (a) (upper panel) C2H abundance in Model CO-1 
compared with our standard model as a function of Ay at four 
times: 10^ yr (dotted lines), 10^ yr (dot-dashed lines), 5 X 10^ yr 
(dashed-lines) and lO'^ yr (solid lines); (b) (lower panel) as in a) 
but for Model CO-2. 



model are relatively close, generally less than an order of 
magnitude, to those of our standard model, the early- and 
late-time behaviour of most of the species is opposite to that 
of the first model, i.e. an even closer agreement to our stan- 
dard model is found at early times. At later times, the lower 
general abundance of elemental C produces smaller abun- 
dances compared to the standard model for all the carbon- 
chain molecules. This effect is specially apparent in the inner 
layers, A„ > 3, with factors that range from 2 to more than 
10. Also at later times, the differences in the outer layers are 
usually less than 50 per cent, and in some cases we can see 
a slight increase in abundance, as is the case for C3H2 or 
CeH, but never affecting the abundance peaks of the stan- 
dard model. This slight increase can also be seen for HNC, 
HC3N, and HCO'*'. The abundance of NH3 is larger by fac- 
tors of 1.5-4 at all layers. 



5 COMPARISON WITH OBSERVATIONS OF 
THE HORSEHEAD NEBULA 

One of the main difficulties in the modelling of PDRs is to 
compare the results with actual observations, because it is 
usually extremely difficult to determine the distribution of 
the gas density, and the way the far-UV radiation penetrates 
the cloud. In recent times, the Horsehead nebula has been 
used as a benchmark for PDR models, because it appears 
to be viewed edge-on and thus offers a good opportunity to 
study the structure of a PDR. 

The Horsehead nebula (B33) is located at the west- 
ern edge of the molecular cloud L1630, at the near side of 
the HII region IC 434. The Horsehead is thought to be an 
early Bok globule; i.e., a higher density cloud core which 
is emerging from its parental and more tenuous cloud as 
the UV radiation from the 09.5V s tar a Ori further to 
the w est erodes the surrounding gas (|Reipurth fc Bouchetl 
Il984h . The Horsehead nebula has been extensively stud- 
ied at many differen t wavelengths, from the early maps in 
CS, CO, and [CII] ('Lada, Ballv & Stark' 'l99ll: IZhou et all 
I1993I : iKramer. Stutzk i & Winncwisscr 199(4) to more recent 
observations ranging from the mid-IR to the millimetre 
llAbergel eral] |2002| . l2003l: iPound Reitaurth fc Ball\] l2003l : 
iTevssier et aiTbOO^ : iHabart et all l2005l : iPetv et all boosf ). 
The gas density derived from these observations, either from 
the excitation of molecular lines or from the penetration 

J 



depth of the UV radiation, is 1-4 x 10 cm" "MZhou et al 



ll993l:)Kramer et al.lll996l : lAbergel et al]|2003l ') . iHabart et al 



( 20051 ) found that the gas density follows a steep gradient 



at the cloud edge, rising to wh = 10* cm ^ in less than 0.02 
pc, and at the same time th e gas temperat u re ch anges from 
~ 200-300 K to ~ 10-20 K. lAbergel et al] (|2003l ) derived a 
kinetic temperature from single-dish observations of 30-40 
K, which would correspond to the bulk of the cloud, while 
the w arm UV-illumina ted edge would have higher tempera- 
tures (jPetv et al.ll2005l l. The far-UV incident radiation field 
can be considered to have a mild intensity, x ~ 60 in Draine 
units , and h as also been estimated in the y ^ 20-200 rang e 
( Zh ou et all [l99 3 ; Kra mer et al.l Il996l : lAbergel et aD l2003l ) . 

iTevssier et al] (|2004l ) and lPetv et al] Hooi) studied the 
physical and chemical structure of B33 at millimetre wave- 
lengths, using moderate-angular-resolution single-dish and 
high- a ngular-resolution inte rferometric observations, respec- 
tively. E^ssiereEaD (|2004h observed the emission of a vari- 
ety of molecules (C^^O, C2H, C3H, C4H, C3H2, C3H4, CeH, 
HN^'^C, HC3N, and CS) at three positions in the Horsehead 
nebula, one coinciding with the peak in IR emission (IR 
peak), another, very close by, at the position at the hydro- 
carbon peak emission, and the last one at the position of the 
CO peak. These and previous observations showed a relative 
stratification of H2 and small hydrocarbons: the boundary of 
the H2 emission corresponds to the edge of the PDR, while 
the hydrocarbon peaks, which are well correlated with one 
another, are located farther into the cloud. In addition, the 
CO p eak occurs at a plac e tangential to the hydrocarbon 
peak. ITevssier et al.l (j2004l ) used steady-state PDR models 
with the same physical conditions as ours to explain the 
original observations, with reaction rates from two diff erent 
networks, UMIST95 fMillar Farauhar fc Willacvlll99'ii ) and 
osu. 2003 (jSmith et al.ll2004h . They reproduced some of the 
observed fractional abundances (C^*0 and C2H) and some 
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column density ratios, but their calculated peak fractional 
abundances for C-C3H2 and C4H are low by an order of mag- 
nitude. The discrepancies were attributed to the return of 
small carbon molecules to the gas phase through the sput- 
tering of small carbon grains and PAHs. 

The high-angular-resolution observations of I Pet v et alj 

determined the column density and fractional abun- 
dance of C^*0, C2H, C4H, and C-C3H2 at three positions, 
named 'IR edge', 'IR peak', and 'cloud', in order of their 
depth inside the cloud. The TR-edge' position corresponds 
to the edge of the PDR, where the UV radiation is most 
intense; the 'IR peak' position is located 5 arcsec east, 
where both t h e IR and hydrocarbon emission peaks of 
iTevssier et all (|2004l ) were found; and the 'cloud' position 
is ~ 15 arcsec deeper into the cloud. The gas tempera- 
tu re estim ates for the three positions are 100, 60, and 40 
K. IPetv e t al. ( 2005) also used different models to try to fit 
the more spatially detailed observations but t he main results 
are no t significan tly different from the ones of pRv ssier et al.l 
120041 '). Recentlv. lPetv et all (|2007h published more data for 
CS, HCS^, and 1100"*" at three positions, one approximately 
at the 'cloud' position and two deeper into the molecular 
cloud. 

The results of our standard model whe n steady-state 
is rea ched are genera l ly sim ilar to those of ITevssier et al.l 
l|2004[l and lPetv et al.l (|2005l ) (with either reaction network) 
concerning the details of the PDR zones and the posi- 
tions and heights of the hydrocarbon peaks. We compare 
in Figs [11] and [12] the abundance profiles of our standard 
model at three different times with th e fractional abun - 
dan ces of some n i olecul es observed by iPetv et ah! (120051 ) 
andpRv ssier et al] (|2004l ). taking into account the errors as- 
sociated with the observations. In these figures, fractional 
abundance is plotted against visual extinction, determined 
from position for the observational values by assuming a 
constant density of 2 x 10* cm~^. The peaking of fractional 
abundances towards the outer layers of the PDR in both 
observation and theory is striking. 

Let us first consider the results at steady-state (5 x 10^ 
yr), focu sing initially on th e molecules studied at high reso- 
lution bv lPetv et all (|2005h . The peak fractional abundances 
of C2II, C4II, and C3 H2 in our standard mo del are similar to 
those of the model of ITevssier et al.l (|2004l ). Compared with 
observation (see Fig. 1111) . we find reasonably good agree- 
ment at all three observational positions for C2H (less so 
for the 'cloud' position, although a slight displacement of 
AAv ~ 0.25 would fit the three points at the same time). 
For C4H and C3H2, the value observed at the 'cloud' posi- 
tion is the only one matched by our model, within a factor of 
2 for C4II. The underproduction for the other two positions 
is ~ 10^ - 10=* times for C4H and <50 times for C3H2. 

Fig. | 12l s how the additional s pecies detected by 



ITevssier et all (120041 ) and IPetv et al.l (|2007l ) compared to 
the results of our standard model. We note that their so- 
called IR peak and hydrocarbon peak positions are assumed 
to both be represented by visual extinctions from to 2. To 
convert the actual observation of IIN^"^C and H^^CO^ into 
the normal isotopologue , we assumed a ^^C/^'^C ratio of 62 
(|Langer fc Penziaj|l993l ). As can be seen, the steady-state 
results of our standard model range from good to poor. CS 
and HCS^ are reproduced within a factor of ~ 3, which for 
the first could probably be attributed to the effects of sulfur 




Figure 11. Comparison of the abundance profiles of our standard 
model for C2H, C4H, and C 3H2 at three times with the values ob- 
served bv lPetv et al.l ||2005| ') at the three positions, labelled from 
left to right as 'IR-edge', 'IR-peak', and 'cloud'. The boxes indi- 
cate the uncertainties in the observed fractional abundances and 
in the position (beamwidth). 




Figure 12. Similar to Fig. \TT\ but for C3H, CeH, HNC, and 
H C3N in the 'IR peak' o r 'hydrocarbon peak' positions observed 
by ITevssier et al.l ||2004) and for CS, HCO+, and HCS+ taken 
from lPetv et al.l l|2007fr 



depletion. C3H and CeH are underproduced in the model, 
similarly to the other carbon chains, by factors of ~ 3 and 
~ 6, respectively. The abundances of HCO^ and HNC are 
almost a factor of 10 below the observed value, although 
there are values in the model at high Ay that match the ob- 
servations. We find the largest differences for HC3N, which 
is heavily underproduced in the standard model at steady 
state. 

If we consider earlier times, we find that for times of 
10'* yr and later, most of the time dependence occurs in 
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inner portions of the PDR, such as the 'cloud' position of 
iPetv et al.l (|2005l ) and even larger values of visual extinction. 
For the 'cloud' position, we find that the abundance profiles 
of C2H may be some what better fitt e d at e arlier times. In 
general, the results of IXevssier et al.] l|2004l ) do not extend 
far enough into the cloud to be sensitive to time. There is a 
large time dependence for the abundance of HC3N but this 
occurs at high values of visual extinction for the most part. 
Since variations in the CO initial abundances discussed in 
Sect. l4.3l mainlv affect abundances in inner layers, with little 
or no impact on the outer peak abundances, these variations 
do not greatly affect the level of agreement between model 
and observation. 



5.1 Other physical parameters 

We used the alternative models discussed in Sect. |4]to de- 
termine if there is a set of physical parameters that agrees 
better than the standard model with the abundances de- 
rived from the observations in the Horsehead nebula. We 
also ran additional models that looked at intermediate val- 
ues of the parameters or that changed initial conditions. 
We did not find a unique best-fit model. Since differe nt au- 
thors (e.g., iPetv et al.1 l2005l : iGoicoechea et al.l l2007h have 
proposed different elemental abundances to try to explain 
the abundances in PDRs, we have also changed the relative 
elemental abundance, mainly by increasing the abundance 
of C, but also of O, S, and N. An initial carbon abundance 
twice of that in our standard model shows some of the bet- 
ter agreements for C4lf (for the innermost position), CeH, 
and especially CS and IfC S^. 'Carbon-r i ch' ab undances sim- 
ilar to the ones used by I Millar et al.l l|200(]| ) overproduce 
the abundances of most of the carbon-containing molecules. 
Different sulfur abundances can be used to fine tune the 
amounts of CS and HCS^. We also explored the effects of 
different values of the cosmic ray ionization rates. We used 
values of 0.5, 10, 50, and 100 times the standard value of 
Typically, higher values of the cosmic-ray photodissocia- 
tion rate (> 10 — 50) show better agreements for many of 
the observed molecules, especially combined with a slightly 
higher density and temperature and/or slightly less intense 
radiation field. 

5.2 Age of the PDR 

The time dependence of our models, as is also the cas e 
for the time-dependent model of UCL (|Bell et al.l l2005h . 
can in principle be used to try to set an age for a dense 
PD R, given a favourabl e geometry. For the Horsehead neb- 
ula, iPound et all l|2003l ) argue that, from its size and veloc- 
ity gradient, the formation time scale could be as short as 
a few X 10^ yr, and the object is likely to be destroyed in 
5 X lO'' yr. These time scales are shorter or comparable to the 
time needed by our standard model to reach steady-state. 
But, in trying to distinguish between steady-state and earlier 
times by chemical abundance profiles, we must keep in mind 
some formidable obstacles. First, the uncertainties in the 
determination of the PDR geometry and density structure 
make such an analysis very difficult because they obscure 
purely temporal effects. Secondly, there is a dependence of 
early-time results on what is chosen for initial abundances. 



Thirdly, the abundance profiles of the many hydrocarbons 
that peak in the outer layers of the cloud show greater sensi- 
tivity to time in the inner layers of the PDR because steady- 
state takes more than 10 times longer to reach in the inner 
slabs. Unless the PDR is very young, these inner slab abun- 
dances are likely to be much smaller than peak values and 
thus hard to observe. Moreover, most of the more abundant 
molecules that peak in the deep slabs of the cloud also reach 
values close to those of steady-state in a relatively short time 
and may not help much in the determination of age. Some 
molecules, such as HC3N, suffer large variations in fractional 
abundance in the same span of time and in the position of 
the peak abundance, and these could be useful to carry out 
age estimates. In some of the models we ran, we find better 
agreements for some molecules at times 0.4 — 3 x 10"" yrs, 
but given the uncertainties in the observations and in the 
density structure of the cloud, we need more observations 
of different positions and molecular transitions to have a 
clearer view of the structure of the PDR and to be able to 
make a more precise estimate of its age. In particular, we 
would need observations of the layers with Ay > 4, where 
wc find, in our model, the largest variations of abundance 
for most of the observed molecules. 



6 SUMMARY 

We present the results of a time-dependent model of a dense 
photon-dominated region created after a star has formed in- 
side a dark cloud. We have assumed for our standard model 
an homogeneous cloud with density n = 2 x 10* cm~^ and 
gas temperature T = 40 K, on one side of which impinges 
an FUV radiation field x = 100 times the ISRF. 

We have found that the more unsaturated and/or sim- 
pler hydrocarbons, in particular C2H, C4H, or CeH, have 
peak fractional abundances in the outer layers of the PDR, 
at Av ~l-2, at times as early as a few xlO* yr, while 
more saturated hydrocarbons, such as CH4, have their peaks 
deeper into the cloud at these times, around Av ~ 5. The 
outer and intermediate layers reach steady-state rapidly, at 
times shorter than a few x 10^ yr, while for the inner layers 
it takes almost an order of magnitude longer. 

The molecule CS and the cyanopolyynes also peak 
nearer to the surface of the cloud, at Av 1.7-2, at times 
following a few xlO'' yr, while HCN, HNC, HCO+, N2H+ 
and NH3 peak deeper into the clouds, from Ay = 5.5 on- 
wards. This result is very interesting because it reproduces 
in space a similar distribution of molecules to the one found 
in time for a homogeneous gas-phase model. The external 
radiation field seems to keep the chemistry 'young' in the 
outer layers of the cloud, while only the more shielded slabs 
develop a 'late-time' chemistry. 

We have run several models with modifications of some 
of the parameters of the standard model, varying the inten- 
sity of the radiation field, the temperature and density of the 
cloud, and the initial distribution of fractional abundances 
in the gas. For the three models run with non-standard scal- 
ing factors of the radiation field of 10, lO'', and 10*, the ef- 
fects on the peak fractional abundances of hydrocarbons are 
not large, typically about an order-of-magnitude. At large 
values of the radiation field, all molecules are swept out of 
the outer regions, while only the gas at Ay > 9 — 10 is 
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completely shielded from variations in the radiation field. 
One salient result of increasing the density is that for hy- 
drocarbons that peak in the outer layers, the abundance 
profile when plotted against distance has a smaller region of 
large fractional abundance. Another salient feature is that 
the peak fractional abundances of selected hydrocarbons and 
cyanopolyynes increase with density in the range 10'' - 10® 
cm~^ while other species show little change or even a de- 
crease. Variations in peak fractional abundance as tempera- 
ture is changed tend to be less significant than with density 
for temperatures in the range 10-200 K. We also explored 
the effects that different initial fractional abundances have 
on the results of the model by running a model that starts 
with atoms in the gas, finding that the evolution of frac- 
tional abundances in time is very different from the standard 
model, although the end result (steady state) is, of course, 
the same. 

In a comparison of the results of our standard model 
with recent high-angular resolution observations of the 
Horsehead nebula, we were able to reproduce, within the 
estimated observational uncertainties, the peak fractional 
abundance and spatial dependence of C2H; not so with C4H 
and C3H2 which show much worse agreements at the out- 
ermost positions. For the rest of the molecules, studied at 
moderate- and high-angular resolution, the agreement to the 
observations go from good (for CS and HCS^) to poor (es- 
pecially for HC3N). We have also used additional models 
to explore how different physical conditions can improve the 
agreement with observations. Even if some models did repro- 
duce better some of the molecular abundances, we could not 
find a single model that was able to fit all the observations 
at the same time. 

We also ran two models in which the initial CO abun- 
dance was changed from its 'standard'value of 7.3 x 10~^ 
with respect to total H. First, we reduced this value to that 
measured in TMCl-CP (4x10"^^), while conserving the total 
elemental abundances of C and O by increasing the atomic 
abundances of carbon and oxygen. Secondly, we used the 
reduced CO abundance with no attempt to maintain our 
standard elemental abundances of C and O. In both in- 
stances, little change was seen for carbon-chain species from 
the standard model at times near steady state. At earlier 
times, however, the first modification produces strong en- 
hancements in carbon-chain molecules for the inner layers 
while the second modification produces smaller diminutions 
of these abundances that are even smaller than at steady 
state. The status of the agreement with observations in the 
Horsehead nebula is unaffected by changes in the initial CO 
abundance, since the changes affect the inner layers, which 
have not been observed. 

Recently, there have been some new results that might 
help to solve the disagreement between observation and the- 
ory. First, Rimmer and Herbst (in preparation) have re- 
calculated the cosmic-ray ionization rate as a function of 
depth into clouds. We found in our models that a differ- 
ent value of ( may dramatically change the abundances of 
the species. Secondly, there is evidence that calculations in- 
cluding negative ions may indeed increase the calculated 
abundances of neutral carbon chain species. These calcu- 
lations are currently being done by R. Ni Chiumin, Tom 
Millar, Nana s e Har ada, Eric Herbst and others. Thirdly, 
iRoueff et al.l (|2007l ') have confirmed that the Horsehead 



Nebula is rather clumpy by their analysis of deuterated 
molecules. We intend to use these new calculations in further 
work to find out if they are able to reproduce the abundances 
of the molecules observed in the Horsehead nebula. 

From the currently available observational data, it is 
difficult to obtain a definitive age of the Horsehead nebula. 
We would need more observations, especially of layers deep 
into the cloud (Av-^4), where the abundances of species 
such as hydrocarbons and cyanopolyynes are more sensitive 
to changes in time. 
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APPENDIX A: CALCULATIONS MADE WITH 
THE NEWEST RATEFILE OSU_03_2008 

The calculations discussed in the body of the paper were 
performed with the osu.2003 network for a variety of rea- 
sons. We have also used the latest version of this network 
without molecular carbon chain anions (osu_03_2008) to de- 
termine if there are any substantial changes in the standard 
version of our model. Some salient differences found are as 
follows: 

(i) The abundance of molecular hydrogen is unchanged, 
and there is very little change in the abundance of atomic 
hydrogen ( ~ -1-3 per cent at Ay ~ 2, ~ —3 per cent in the 
innermost layers). The crossing of the H and H2 abundances 
is at Av ~ 0.04. CO abundances are ~ —5 per cent up to 
Av ~ 2, while C shows the opposite behaviour. Abundances 
of C"*" at layers Ay > 4 are a factor of ^ 25 per cent less. 

(ii) Differences between abundances arc quite small for 
most molecules (from a few percent to a factor of 2) at steady 
state, especially for the more abundant hydrocarbons, for 
which they are typically much less than a factor of two, 
except at middle layers {Av ~4-7), where these differences 
can go from 25 per cent less for C2H to a factor of 2 less 
for C3H2. The main apparent effect is a slightly lower inner 
abundance peak for some of these molecules. Other neutrals 
show very small differences, such as HON (abundances 20 
per cent larger or lower overall, except a factor of 2 less at 
Av ~ 4) or H2CO (with almost no changes, less than 5 per 
cent). HCO"'" shows up and down variations of about 20 per 
cent at steady state. Some less abundant hydrocarbons and 
other low abundance species may show larger differences at 
some positions but always at positions where they have very 
low abundance. 

(iii) N2H"'" is generally larger with the newer network, 
with differences up to a factor of five at some outer lay- 
ers at steady state, whore the ion is loss abundant. At the 
inner layers, the differences range from 50 per cent to 20 per 
cent from early times to steady state. Ammonia does not 
show this behavior; its abundance is at most 70 per cent 
larger at middle layers at steady state, but small differences, 
less than 10 per cent in the rest between the two networks 
at steady state (differences of ~ 40 per cent at early times). 

(iv) CS and HCS^ have lower abundances with the newer 
network, with maximum differences a factor of 3-5 for CS 
and 2-15 for HCS"*", especially for the outer layers, and for 
t > 5 X 10^ yr. At the same time, SO and SO2 show larger 
abundances for layers Ay > 6, of ~ 70 per cent and a fac- 
tor of about 4, respectively, at steady state (with smaller 
differences at earlier times). 



(v) HC3N has a larger abundance at steady state, by al- 
most an order of magnitude, for the outer layers, where it 
is less abundant, and of ~ 2 in the inner layers. At earlier 
times, the abundances at the inner layers can be a factor of 
~ 5 larger. The other cyanopolyyncs tend to have slightly 
lower abundances overall, except a drop at Ay 5 that can 
be about an order of magnitude, but they have very low 
abundances in both cases. 

(vi) In comparison with the observational results for the 
Horsehead Nebula, there arc few significant differences. CS 
and HCS"*" would be the most affected molecules, specially 
the latter. Nonetheless, the HCS"*" abundances at Av^ ~ 2-4 
would still be inside the error bars of the abundance deter- 
mined by the observations. 



